Sharing polarization within quantum subspaces 
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Given an ensemble of n spins, at least some of which are partially polarized, we investigate the 
sharing of this polarization within a subspace of k spins. We assume that the sharing results in a 
pseudopure state, characterized by a single purity parameter which we call the bias. As a concrete 
example we consider ensembles of spin- 1/2 nuclei in liquid-state nuclear magnetic resonance (NMR) 
systems. The shared bias levels are compared with some current entanglement bounds to determine 
whether the reduced subspaces can give rise to entangled states. 
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I. INTRODUCTION 

Liquid-state nuclear magnetic resonance (NMR) tech- 
niques have proved a convenient test-bed for implement- 
ing quantum algorithms 0, 0, 01 El IS ■ However, a con- 
ventional NMR quantum computer starts off in an initial 
state which is very similar to a maximally mixed state, 
due to the tiny population differences between the energy 
levels of the spin- 1/2 nuclei which act as the qubits 0>l3- 

Various approaches to this initialization problem have 
been suggested. Firstly, it is possible to convert thermal 
states into pseudopure states Q, Q: this approach has 
been adopted in the vast majority of NMR implemen- 
tations of QIP, but is not scalable J8|. Secondly, non- 
Boltzmann distributions can be prepared, for example 
with para- hydrogen 0,113,113, 

giving almost pure states 
which lie above the entanglement threshold lia| . Thirdly, 
computational schemes have been described |l3l IbH ITq 
which concentrate the small polarization available on a 
large number of spins into a smaller subspace. 

Here we investigate the effects of polarization sharing 
on achieving states of useful purity. The schemes men- 
tioned above concentrate polarization onto a smaller sub- 
set of spins, but we study the effects of spreading the po- 
larization over a larger subset of spins. We assume that 
the polarization sharing results in a pseudopure state 
with a purity characterized by a bias parameter 8 that 
we define below. The bias is then compared with known 
entanglement bounds derived for pseudopure states. 



II. PSEUDOPURE STATES AND 
ENTANGLEMENT 



In the high temperature approximation, the initial 
state of a (homonuclear) NMR quantum computer is 




(1) 



where n is the total number of spins, l n is the identity 
matrix of order 2", B = —tuo/kBT is a Boltzmann factor 
and o~j Z is the z Pauli matrix for spin j. A standard 
implementation of an algorithm requires the quantum 
computer to start off in a pure state, characterised by 
having a single non-zero eigenvalue, of size one. Clearly, 
the thermal state Q has many different eigenvalues and 
so to prepare pure from thermal states, a non- unitary 
process must be employed at some stage. Cooling to the 
ground state is an example of such a process, but is in- 
effective in NMR due to the small energy gaps involved. 
Pseudopure fc-qubit states have an eigenvalue spectrum 
between these extremes, having one large eigenvalue and 
2 k — 1 degenerate smaller eigenvalues. This means that 
thermal states can be converted into pseudopure states 
by particularly simple non-unitary processes, such as av- 
eraging 2 fe — 1 population terms. 

A general fc-qubit pseudopure state x takes the form 



(2) 
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with dynamics identical to those of the corresponding 
pure state \yb){ip\. (We explicitly use the symbol x to 
denote pseudopure states.) The state is characterised by 
a single parameter, 6, which we call the bias of the state. 
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Another useful description is to write out the pseudopure 
state explicitly in its eigenbasis with the eigenvalues in 
descending order to give 



(2-1) terms 



Xk,8 = Dg[/, 



1-/ I"/ 



2 k _ i ' 2 k - 1 ; 



2 fe -l J 



(3) 



where Dg[. . .] denotes a diagonal matrix and / is both 
the fractional population of the desired state \4>) and 
the largest eigenvalue and is related to the bias by 



5 = 



2 k f-l 
2 k - 1 ' 



(4) 



A pseudopure state will always be separable if S is less 
than some critical value, and Braunstein et al. |l6j | have 
described upper and lower bounds on separability. A 
state of the form given in Eq.|2]was shown to be explicitly 
separable for sufficiently small biases such that 



5<8t 



1 



1 + 2 2fc -! 



(5) 



and we say that these states belong to the region S. If 
the bias exceeds their upper bound 



8 > 5 U 



1 



1 + 2* 



(6) 



the state is said to lie in the entangled region E, or more 
appropriately the entanglable [T^ region. In the region in 
between, ES, it is not known whether entangled states 
can be prepared. These bounds were subsequently im- 
proved, shrinking the size of ES. For example, Gurvits 
and Barnum |18| have tightened the lower bound to 



8, = 



2(6) fc / 2 



(7) 



Current NMR implementations with Boltzmann initial- 
ization use states lying in the separable region S and 
it seems possible to enter the region ES only by em- 
ploying more qubits. For example, using the Gurvits- 
Barnum lower bound with a typical Boltzmann factor of 
B = 10 -5 , the state will cross over into ES for k > 41. 

Similar themes have been taken up by other authors 
as well. Yu et al. have shown |17| that by using unitarily 
transformed thermal states in place of pseudopure states, 
the ES to E transition can take place with a smaller num- 
ber of qubits: making pseudopure states from thermal 
states involves convex mixing and, therefore, decreases 
the likelihood of entanglement. Other researchers have 
shown that entanglement can also exist in a 2 x N di- 
mensional quantum system, when only the qubit is pure 
and the N dimensional system is in a highly mixed state, 
such as a two level atom interacting with a high temper- 
ature field. 

Here we investigate the sharing of polarization within 
a pseudopure quantum subspace. We consider an initial 



state 



Pn,p,i 




(8) 



where n is the total number of qubits of which p qubits 
are in the state, 



(l + cr)/2 

(l-cr)/2 



< a < 1, 



(9) 



having a polarization a and the remaining n—p qubits are 
maximally mixed. Our goal is to compute the achievable 
bias in a k qubit pseudopure state where k < n, i.e., 
we are interested in the transformation p n ,p,a l— > Xk,s- 
We achieve this in two steps: (a) convert p n ,p,a into a k 
qubit state using a partial trace operation, and then 
(b) convert Tfc into a pseudopure state \k.s using cyclic 
averaging. 

For simplicity consider a state p n ,p,i with all the p 
qubits being perfectly polarized, 



Pn,p,l = | {&>|0)(0| ® I U/2 



(10) 



= -L_D g ^T7..., i,o,o,. ..,o]. 



2 n—p 



and consider the case k = p. Partially tracing |2jj out 
the n — k qubits from p n ,p,i leaves us with the k qubit 
reduced state 



2-1 terms 



r fe = Dg[l,0,0,.7T0]. 



(11) 



This is also a pseudopure state Xk,i with maximum 
achievable bias of one. The partial trace operation re- 
duces the dimensionality of the system and for a diagonal 
state, is equivalent to taking partial sums of consecutive 
eigenvalues along the ordered diagonal. Physically, it just 
corresponds to ignoring the n — k qubits. In the case of 
k = p, the reduced subspace after the partial trace opera- 
tion is already a (pseudo)pure state, with 8 — 1. This will 
also be the case when k < p. However, for p < k < n, the 
reduced state will have several different eigenvalues. 

We can explicitly generate the desired pseudopure 
state by writing the mixed state in its ordered eigen- 
basis and then averaging over cyclic permutations of the 
2 k — 1 trailing elements. This procedure (which can be 
considered as a generalisation of the twirl operation plf) 
is experimentally implementable by exhaustive temporal 
averagin g |22l|: more efficient procedures are also avail- 
able [22L |2j| . A non- unitary process cannot increase the 
maximum eigenvalue of a state r, but by definition, a 
pseudopure state has only one non-degenerate eigenvalue 
which must coincide with the maximum eigenvalue of r. 
Cyclic averaging leaves this eigenvalue unchanged and 
therefore extracts the maximum bias 8. 
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FIG. 1: The effect of sharing the polarization of p pure [(a)- 
(e)] or partially pure [(f)-(i)] spins on quantum subspaces. 
Starting with p pure spins (a), we can pick k < p pure spins 
(b) and extract all the available polarization. The polariza- 
tion can also be spread onto a bigger subspace while retaining 
the entanglability (c). Sharing the polarization within even 
higher dimensional subspaces decreases the resulting bias and 
therefore the possibility of entanglement, (d) and (e). As the 
size of the extracted subspace increases, the bias drops expo- 
nentially. Smaller biases are shown by lighter tones of grey. 
For p partially polarized spins (f), we can concentrate the bias 
onto smaller subspaces (g) or spread it onto bigger subspaces, 
(h) and (i), with a corresponding decrease in the bias reducing 
the likelihood of entanglement. 



Algebraic manipulation shows that the maximum 
eigenvalue of Xk,5 is given by the formula 



/ = 



Tl'(pri,p.<rPri,fc,l) 

Tr(< M ) 



1 



zn:T^{Pn, P ,<Tpn,k,i) (12) 



which is the overlap between the states (jHJl and (|10J) . 
Note that we do not consider here precisely how the po- 
larization sharing procedures might be implemented in 
practice, but simply determine the limits. Our results 
are summarised in Fig.n an d exemplified below. 



III. PERFECTLY POLARIZED QUBITS 

Considering the form of Equation 1101 (n spins out of 
which p are perfectly polarized and the rest are maxi- 
mally mixed), we see that / is given by 




for k > p 
for k < p 



and from Equation 01 the achievable bias is 

for k > p 
for k < p. 




(13) 



(14) 



For k < p, all the purity can be extracted as this is 
equivalent to picking out the pure spins from a set of 



0.01 




FIG. 2: (Color online) Sharing the polarization of 2 and 4 pure 
spins. The upper shaded region is E and the lower shaded 
region is S, whereas ES is in between. The border between 
two regions belongs to the less entangled region. The regions 
are intended to guide the eye, and are only shown for 2 or 
more spins, where the concept of entanglement is valid. 



pure and maximally mixed spins. The case k > p is more 
interesting as this involves distributing the polarization 
of p spins over a larger subspace and is the polarization 
sharing that we shall be mainly interested in. Figure |21 
shows the results for p pure spins, with p = 2 and 4. For 
k > p, the bias falls off exponentially as k is increased. 

Examining Equation 1101 also shows that adding extra 
qubits (increasing n) neither hinders nor helps the extrac- 
tion of purity onto qubits; it seems that the only role of 
n is that it limits the number of spins in the pseudopure 
subspace, as evidently k < n. Of course these additional 
qubits may assist the implementation of the manipula- 
tions necessary to prepare the pseudopure state, but they 
have no ultimate effect on the extracted bias. 

The biases of the fc-qubit pseudopure subspaces can 
now be compared with the entanglement bounds, but for 
the possibility of entanglement to arise in the first place, 
the pure component in the pseudopure state must be en- 
tangled. For example, the direct product 2-qubit state 
00} {00 1 is pure, but is clearly separable. However, one 
can always find a unitary operator that converts an ar- 
bitrary pure state into any desired pure target state, be 
it entangled or otherwise. Any pure state is therefore 
unitarily equivalent to an entangled state of the same 
dimensions, and we can assume that our fc-qubit pseudo- 
pure state comprises a pure, entangled component mixed 
in with the maximally mixed state. This is the justifica- 
tion for using the term entanglable for the region E, as 
mentioned in the previous Section. 

For k < p, the fc-qubit subspace is obviously in the 
region E but for k > p, increasing k leads to an expo- 
nential drop in the extractable bias 5. At some critical 
k we step into the region ES. This transition from E to 
ES takes place where 5 = 5 U , which can be identified in 
Figure |21 as the point at which the bias curve crosses into 
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the region ES, and occurs at k c = 2p. What this means 
is that if we start with p perfectly polarized spins, then 
for k < 2p our extracted system will be in E. However, 
for higher values of k we shall be in ES. We also note 
that as a result of this polarization sharing we will never 
enter S; this is clear from the slopes of the lines in the 
figures, and a simple proof is given in the Appendix. 



IV. PARTIALLY POLARIZED QUBITS 

We now generalize the above example to consider p 
spins that are not perfectly polarized, each having a uni- 
form non-zero polarization < a < 1. The n-qubit state 
Pn, P ,a- is now given by 



Pn.p.a 




j=p+l 
(2"- p p Ci) terms 




where l Cj — l\/ — j)\ is the binomial coefficient. If 
k > p, calculating / involves summing only the first 2™ 
terms, each being of the same size ((1 + a)/2) p . This is 
best illustrated with a numerical example. Consider a 
state with n = 4 and p = 2 



4 terms 



4 terms 



Pi,2,a = ^ Dg[ 

(16) 

and suppose we want to share the purity over a 3-qubit 
subspace (k = 3). Calculating / will only involve a par- 
tial sum of the first 2 n ~ k = 2 terms in (|16|l . each of these 
terms being (1/4)((1 + er)/2) 2 . Similarly if k = 2, we 
need a partial sum over the first 4 terms and for k = 4, 
we need to consider just the first term. In each of these 
cases only the leading terms are involved in the partial 
sums, and we can derive a formula for / 



f = 2- k (\ + ay 
with a corresponding bias 

(l + of-l 



6 = 



2 k - 1 



(17) 



(18) 



If k < p we are concentrating polarization onto a sub- 
space smaller than the original, and can expect to extract 
a higher bias, in this case the overlap and bias must be 




FIG. 3: (Color online) Sharing the polarization of impure 
spins. The regions E, ES and S are highlighted in the same 
way as in Fig. The figure demonstrates polarization shar- 
ing for 5 pure and 5 impure spins. The effect of sharing from a 
single weakly polarized spin is also shown; when shared over 6 
or more spins the state moves from being provably separable 
into the region where it is unknown whether or not it may 
be entangled. Note, however, that states currently placed in 
ES may subsequently move to S if lower bounds are further 
tightened. 



calculated explicitly on a case by case basis. We do not 
consider this polarization concentration further. 

As with perfectly polarized spins, sharing from im- 
pure spins also entails an exponential decrease in the ex- 
tractable bias as the size of the subspace k increases. The 
critical size k c for impure states at which the state is no 
longer provably entangled is now given by 



fc c = \2p 



ln(l + a) . 
ln~2 



(19) 



where \x~\ represents the next higher integer to x (only 
integral numbers of spins are possible). It is straightfor- 
ward to see that with p partially polarized spins, polar- 
ization can be shared over a smaller number of qubits 
than when the p spins are ideally polarized, as shown in 
Figure |21 



V. CONCLUSIONS 

In this paper, we have shown that given p perfectly 
polarized spins, we can share polarization among 2p — 1 
spins, such that their bias keeps them in the provably 
entangled region. For example, for p = 2 pure spins, we 
can spread the polarization over 3 spins, while "protect- 
ing" the entanglement. Similar results can be achieved 
when the initial spins are not completely polarized: from 
we find that for a > (exp(31n2/4) - 1) » 0.682, 
it is possible to spread the polarization of 2 qubits onto 
a higher spin subspace. Our results suggest that polar- 
ization sharing is of limited value with only a single pair 
of protons from para-hydrogen added onto our substrate 
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molecule |9( , but show more promise for a higher number 
of molecules M. Preparing such molecules is in principle 
possible and we are currently investigating approaches 
for such multi-qubit systems. 
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APPENDIX 

When sharing polarization from p pure qubits, the sys- 
tem always lies within the regions E and ES, never en- 
tering the explicitly separable region S. To prove this we 



must show that the extracted purity is always greater 
than the Gurvits-Barnum bound, 



2 P - 1 



> 



2 k - 1 2(6) fc / 2 ' 



V k > p. 



(A.l) 



Proof. The left hand side of Eauation lA.il will be a mini- 
mum when p = 1 , and so it suffices to prove this extreme 
case: 



1 



> 



2 k - 1 2(6) fe / 



(A.2) 



Now the L.H.S. of (|A.2|) is clearly greater than l/2 fe . So 
the inequality will be true when the ratio of l/2 fc and 
right hand side of (|A.2(1 is greater than unity. This ratio 
(2/3) 1_fc / 2 is greater than one for k > 2, whereas l|A.2() 
can also be shown to be correct for k — 2 by explicit 
calculation. □ 
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